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£f~) ' We study the formation of fractal structure in one-dimensional many-body systems 



with attractive power-law potentials. Numerical analysis shows that the range of the 
index of the power for which fractal structure emerges is limited. Dependence of the 



growth rate on wavenumber and power-index is obtained by linear analysis of the 



collisionless Boltzmann equation, which supports the numerical results. 
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(N ! I. INTRODUCTION 

in 
o 

Formation of spatial structures is an interesting and important phenomenon in nature. 

I : n n 

It is seen over a wide range, from protein folding in biological systems |2j to large-scale 



, ^ ' structure in the universe 

X 
S3 



The theoretical origins of such structures are quite important 
and will be classified into several classes. One of the most interesting areas within the 



field of dynamical systems is that some remarkable structure and organization is created 
dynamically by the mutual interaction among the elements J, ^]. 

Recently we have discovered that spatial structure with fractal distribution emerges spon- 
taneously from uniformly random initial conditions in a one-dimensional self-gravitating 

\ 



system, that is the sheet model 



What is noteworthy in this phenomenon is that the 
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spatial structure is not given at the initial condition, but dynamically created from a state 
without spatial correlation. Succeeding research clarified that the structure is created first 

Q 

in small spatial scale then grows up to large scale through hierarchical clustering p|, and 
the structure is transient 8| . It is quite interesting that some remarkable spatial structures 
are emerged instead of monotonous thermal relaxation in Hamiltonian system. 

The emergence of fractal structure is a typical example that systems of many degrees of 
freedom are self-organized by dynamics themselves. Hence to clarify its dynamical mecha- 
nism is very important subject toward understanding physics of self-organization of matter. 

A way to clarify the dynamical mechanism is to know which class of pair potential can 
form the fractal structure. Here we note an important fact that fractal structure does not 
have characteristic spatial scale, nor does the potential of the sheet model, since the pair 
potential is power of the distance. Hence the scale free property of potential may be a key- 
stone to understand the dynamical mechanism. The question is if the fractal structure can 
be formed in not only the sheet model, but also other systems with power-law potentials. 

In this paper we study the possibility that the fractal structure can be formed in more 
general systems without characteristic spatial scale which is extended from sheet model. 
Here we adopt the model as the system with attractive power-law potentials. At first we 
examine the formation of the fractal structure by numerical simulation for various values of 
the power index of the potential. Next we perform linear analysis to consider the numerical 
results. 

In SecOwe introduce many-body systems with power-law potentials. In Sec JIIII we review 
the formation of power-law correlation in the sheet model. In Sec lIVI we carry out numerical 
simulation. In Sec|V]we analyze linear perturbation of the collision-less Boltzmann equation. 
The final section is devoted to summary and discussions. 
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II. MODEL 



We consider the model where many particles with an uniform mass interact with purely 



attractive pair potential of power-law, which is described by the Hamiltonian 

N 2 N ~ 1 N 



Q 



H = K + U 



Ef+EE 
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(1) 



i=l i=l j>i 

where Xi and Pi are the position and momentum a particle, respectively. The first term 
is kinetic energy and the second term is potential energy. For simplicity in this paper we 
consider the system where motion of particles is bounded to one-dimensional direction. 

For the special case a — 1, the Hamiltonian (0) applies to a system of N infinite parallel 
mass sheets, where each sheet extends over a plane parallel to the yz plane and moves along 
the x axis under the mutual gravitational attraction of all the other sheets. The Hamiltonian 

nnnn 

of the sheet model jlOL II 1L 1121 1131 ] is usually written in the form 

N 2 N ~ l N 



* = Ef + EE 



In the previous letters 



Xi Xj | . (2) 

I i=l j>i 

jj, |8| we investigated the time evolution of the sheet model (J2J) 
to show that fractal structure emerges from non-fractal initial conditions. In the Sees. II VI 
and |3 we will clarify the relation between the power index of the pair potential a and the 
formation of the fractal structure by investigating time evolution of more general model (JT|) 
numerically and analytically. 



III. FORMATION OF FRACTAL STRUCTURE IN THE SHEET MODEL (a = 1) 



Before we investigate the time evolution of the model ([!}, we briefly review our previous 
works f], 2, ^ on the model (J2J). We found that fractal structure is formed in the sheet 
model (J2J from uniformly-random initial conditions f], 7|. In Figure ^ we show a typical 
example of such fractal structure and a snapshot of the /x-space configuration. 



4 



IV. NUMERICAL SIMULATION 

In this section we examine by computer simulations if the fractal structure is formed in 
the systems ((I}. Here we carry out the case a > 1 to avoid divergence when two particles 
come close. 

The non-fractal initial conditions where fractal structure is emerged are characterized 

n 

typically by those of virial ratio 2E kin / E pot = in the sheet model p. (Spatial distribution 
is set to be random.) These state of zero velocity dispersion corresponds to the limiting 
case of zero thermal fluctuation, and is called cold-random condition. Therefore we adopt 
this cold-random initial condition to investigate the time evolution for various value of a. 
In our simulation we use 4-th order of the symplectic integrator with a fixed time step 
At = 27r/10 4 and iV = 65536 particles. In what follows we show numerical results for two 
typical examples: a = 1.125 and a = 1.5. For other values of a behavior of the systems 
varies gradually in accordance with the change of a. 

A. The case a = 1.125 

At first, we consider the case the interaction force deviates slightly from sheet model. 
Here we show the numerical results in the case a = 1.125 in FigsEl and El We display 
particle distribution in (x, v) space (/z-space) in FigEl In the course of time evolution we 
see that many whirlpools nest in the /z-space to form the hierarchical structure. In FigJ3]we 
show a box counting dimension of the spatial distribution. We can see that the dimension 
is D ~ 0.83 (FigOl). Therefore our numerical results suggest that the fractal structure is 



formed. These behaviors are similar qualitatively to the sheet model, a = 1 
that fractal structure can be formed as well as the sheet model. 



. We find 
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B. The case a = 1.5 

Next, we consider the interaction deviate further from sheet model. Time evolution 
change qualitatively as a increases. Here we show the numerical results in the case a = 1.5. 
We display particle distribution in /i-space in FigEl and in FigJ3] we show a box counting 
dimension for spatial distribution. Differently from the case a = 1.125, a single spiral is 
rolled up in /z-space. Therefore fractal structure is not formed. 

We can summarize these numerical results in this section that the fractal structure can 
be constructed as well as the sheet model, but range of a for which the fractal structure is 
created is limited; it can not be constructed for large value of a. 



V. ANALYSIS OF LINEAR PERTURBATION 

In this section we clarify the physical reason analytically why fractal structure can not be 
constructed for the large value of the power index of the potential in numerical simulation in 
the previous section. The formation of the fractal structure occurs at the relative early stage 



in the whole-evolution history Then it is instructive for understanding the mechanism by 
which the structure is formed to know the qualitative properties of the short-term behaviors 
by linear analysis. 

In this section we derive the dispersion relation from the collisionless Boltzmann equation 



14. 



id- 



(CBE), which describes the growing rate of the linear perturbation 

A. Collisionless Boltzmann equation 

CBE is defined by 

di +P dx + J dx'F(x-x')\J dp'f(x',p',t)\ —j f(x,p,t) = 0, 
where F is 2-particle force which is related with a pair-potential U by 



(3) 



^/ x dU , s 

= "TP < 4 > 
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and / is the one-particle distribution function. For simplicity we consider the system is 
extended in infinite region — oo < x < oo. It is clear that the state of uniform spatial 
density with an arbitrary velocity distribution, 

f(x,p,t) = / (p), (5) 
is a stationary state. We impose the following small perturbation over the stationary state 

f(x,p,t)=f (p) + Sf(x,p,t). (6) 
The linearized equation for Sf is 

g- t +p-^Uf(x,p,t) = - J dx'F{x-x') J dp'Sf(x',p',t)—f (p). (7) 

Now we define the Fourier-Laplace transform by that x is Fourier transformed and t is 
Laplace transformed, that is 

roc poo 

6f(k,p,u)= dte~ twt dxe ikx Sf(x,p,t). (8) 

JO J-oo 

Fourier transform is 

/oo 
dxe* kx 5f(x,p,t), (9) 
-oo 

and 



oo 



F(k) = / dxe ikx F(x). (10) 

J — oo 

Then the Fourier-Laplace transformed equation of Eq. ((7j) is 



where Sk{<jj) is 



e k (co)5f k (cu) = L—5f(k,p,0), (11) 

i{— io + kp) 



^) sl+ ',*A> (12) 
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B. Dispersion relation 



From the inverse-Laplace transform of Eq. ljllj) . 



Sf k (t) = / e"*5f h (u>)— (13) 

Z7T 



oo— icr 



" W e !Ui 4lv -krT^^ )^- ( 14 ) 

.-«r efc(w)i(-w + fcp) 2tt 

Now we continuously move the integration contour to upper half of complex u plane while 
avoiding the singular points. Then the contributions from except of the pole can be ne- 
glected, because of the factor exp(iut) (Im(u) > 0). 

Then the growth rate of each mode of the fluctuation is obtained by the solution of the 
equation 

£fcH = 0. (15) 

Eq. (ll5jl is the dispersion relation. If Eq. (|15j) has the solution where the inequality Im(u) < 
is satisfied, the fluctuation is unstable. 

C. Dynamical stability of systems with power- law potentials 

Now we consider the case that the potential is power-law, the pair-potential U is 

U(x)=A\x\ a . (16) 
Assuming the interaction is attractive, A > 0. a = 1 for "sheet model". The Fourier- 

3 



transformed potential is 



U(k) = / dxe lkx A\x\ a (17) 

J — oo 

, [ f air\ T(a + 1) 1 , . 

= 2A -(sin— ) +1 J (a^0,2,4,---,-l,-3,---). (18) 



2 J \k\ a 

For simplicity, we choose the stationary state as 

fo = n 6(p), (19) 
where n is the number density of particles. 
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The dispersion relation is 

e *(u,) = 1 + 2n A | (sin — ) -f^i | - = 0, (20) 
and u which satisfy the dispersion relation is 

- =- 2 "» A {( sin T)V^r (21) 

When u 2 < 0, the system is unstable. Eq. (}2*T]) can be reduced to 



to 



±iJ2n Q A | (sin ^) T(a + 1)} • |A;| (1 - q)/2 . (22) 



From Eq. (J22j) we can classify the evolution of the perturbation into three types; (i) when < 
a < 1, the growing rate increase monotonously for \k\. (ii) when a — 1, the fluctuations for 
all scale grows at same rate, (iii) when 1 < a < 2, the growing rate decrease monotonously 
for \k\. 

The cold-random condition corresponds to the mixed states of fluctuation with all wave- 
length, so called "white noise". The larger the value of a is, the smaller the growth rate 
of the fluctuation with large wavevector is. This is consistent with the numerical results in 
Sec. IIVI that a large whirlpool is formed in ji space in the case for the large value of a. 

VI. CONCLUDING REMARK 

In this paper we have studied structure formation of many-body systems with power-law 
potentials. Firstly we have investigated structure formation in this model by numerical sim- 
ulation. As results, we have found that behaviors of time evolution are different depending 
on the power index of the potential. When a is slightly above 1, fractal structure is formed 
similar to the sheet model (J. On the other hand, when a > 1.5, fractal structure is not 
formed. 

In order to explain these numerical results, we have also analyzed linear perturbation 
of collisionless Boltzmann equation to derive the dispersion relation which represent the 
growth rate of each mode of the fluctuation. As results, we have found that for large value 
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of a the growth rate of small scale is suppressed. In addition we can explain the sheet 
model (a = 1) is marginal in the sense all scale fluctuations grow at same rate. There is 
a slight difference between the initial condition used in the numerical simulation and the 
stationary state employed in the perturbation, that is, we set N particles in finite region for 
simulation, whereas the unperturbed state fo is infinitely spread. Nontheless we think the 
linear perturbation grabs the core of the instability, especially for short time scale. 



ustering. That 
2j. The spatial 



In our simulation the fractal structure is formed through hierarchical c 
is, clusters are created first in small spatial scale, then grow to large scale 
randomness in the initial condition and finiteness of the number of particles N imply that 
relative fluctuation in mass density is large for small spatial scale. This enhancement of initial 
fluctuation in small spatial scale is probably the "seed" of spatial structure formation. 

Fractal properties of the structures can differ from place to place of the system. Most 
fractal structures are not exactly self-similar but can contain various inner structures. Using 
generalized dimension D q one can unveil details of fractal structures. Here we show the 
generalized dimension of the spatial structure formed in the power-law potential model with 
a = 1.125 in FiglHl We find that the generalized dimension is almost constant (D q ~ 0.83) 
for some range of q at t = 5.0 (FigJHl top). Therefore we conclude that the fractal structure 
formed is mono-fractal at t = 5.0 as far as we observed. 

As discussed in our previous paper jsj, however, the fractal structure is a transient state 
and relaxes finally. We find that the multifractal nature emerges due to the relaxation at 
late times (the bottum of Fig. EJ). The boxcounting dimension Dq has relaxed to almost 1, 
while the correlation dimension D 2 has not yet (the bottom of FigEJ). In other words, the 
relaxation of the boxcounting dimension D is much faster than the correlation dimension 
D2 . The detailed analysis will be a future work. 

We noted that scale free property of the potential may be a key-stone to understand the 
dynamical mechanism of the emergence of fractal spatial structure. In this paper we have 
clarified that scale free property does not immediately imply fractal structure. Potentials 
with different power-index makes different temporal behaviors, and the gravitational system 
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the sheet model ) is a special case in all the 1- dimensional systems with scale-free potential. 
As relevant works, structure frormation in other one-dimensional self-gravitating systems 



in an expanding universe has been studied 13, 18]. They c 



aimed that their system is not a 
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Indications of this behavior 



simple fractal or even a regular multi-fractal, but bifractal 
have also been found for the "quintic" model It is a future work to clarify flic relevance 
between the models comprehensibly. In addision, the dependence of the spatial dimension 
on the structure formation will be also interesting and important subject of future. 
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APPENDIX A: PROOF OF EQ.flTT]) 

In this appendix we lead Eq. fjllj) by the Fourier-Laplace transform of Eq.([7j). In terms 
of the transform (jllj) , Eq. (jJJ) become following; The first term of the left hand side of Eq. 
(J7J) is transformed to 



OO POO 







a 



dt 



-i t=oo 



dte~ lu}t dxe lkx —Sf(x,p,t)=e- tu}t 5f(k,p,t) - (-iu) dte- lujt 5f(k,p,t) 



t=o 



OO 



(Al) 

= -6f(k,p,t = 0) + iuTf{k,p,uj), (A2) 

here u is a complex number and Irmo < 0. The second term of the left hand side of Eq. (J7j) 
is transformed to 

/oo o poo 

dxe ikx p—5f{x,p,t) = -ikp / dte~ luJt 6f{k,p,t) = -ik P 8f(k,p,u). 
-oo & x JO 

(A3) 

Then the left hand side of Eq. (J7J) is transformed to 

(Ihs) = + kp)5f(k,p,u) -5f(k,p,0). (A4) 
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On the other hand, the right hand side of Eq. (JJJ) is transformed to 

poo poo f poo poo o ^ 

(rhs)=- dte~ iwt dxe ikx \ dx'F(x-x') dp' 5 f (x' , p' , t) — f (p) \ 

JO J-oo {J-oo J-oo OV ) 



oo poo 



dye^F(y) \ dp'5f(k,p',uj)^(p) 



oo 



F(k) I dp'5f{k,p',u) d -^{p), 



where 



oc 



Sf{x',p',u)= dte- l ^5f(x',p',t). (A5) 
Jo 

Eq. (JJj) is transformed to 

= iF^rfeyf w ' {fy rmM ) + K^TWi^P,o), ( A6) 

where 

poo 

8fk(u) = / dp5f(k,p,u). 



Integrating both side of Eq. (jA6|) by p, we obtain 

1 + r dp F{k } ^(p)) 5f k {u) = ~] 6f(k,p,0). (A7) 
J-oo H - ^ + kp) ®P J + kP) 

Using Eq.flnj), we obtain Eq. ffTTJl. 
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FIG. 1: An example of fractal structure is formed of the model (J2J) from non- fractal structure. The 
number of particles is N = 2 15 . The left figure represents the \i space at t = 9.375. The right 
figure represents the box-counting dimension. 

FIG. 2: The snap-shots of the [i space for a = 1.125. N = 65536. Time are t = 4.4, 5.0, 9.4 from 
the top to the bottom. 



FIG. 3: Box counting dimension D of the spatial distribution for a = 1.125 and t = 5.0. The 
plus symbols represent our data, and solid and dashed lines correspond to D = 0.83 and D = 1, 
respectively. 
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FIG. 4: The snap-shots of the ji space for a = 1.5. N = 65536. Time are t = 11.0, 15.7 and 23.6 
from the top to the bottom. 



FIG. 5: Box counting dimension D of the spatial distribution for a = 1.5 and t = 23.6. 



FIG. 6: The snap-shot of the fi space for a = 1.125 (left). The generalized dimension D q (right). 
Time are t = 5.0, 9.4 and 37.7 from the top to the bottom. 



